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A deep and detailed examination of 29 classical Cepheids with the Spitzer 
Space Telescope has revealed three stars with strong nearby extended emission 
detected in multiple bands which appears to be physically associated with the 
stars. RS Pup was already known to possess extended infrared emission, while 
the extended emission around the other two stars (S Mus and 5 Cep) is newly 
discovered in our observations. Four other stars (GH Lup, £ Car, T Mon and 
q X Cyg) show tentative evidence for extended infrared emission. An unusual elon- 

gated extended object next to SZ Tau appears to be a background or foreground 
ci ', object in a chance alignment with the Cepheid. The inferred mass loss rates 

upper limits for S Mus and 5 Cep are in the range from 10" 9 to 10 -8 M yr _1 , 
with the upper limit for RS Pup as high as 10~ 6 M yr" 1 . Mass loss during post- 
\Q \ main-sequence evolution has been proposed as a resolution to the discrepancy 

^ between pulsational and dynamical masses of Cepheid variable stars: dust in 

the lost material would make itself known by the presence of an infrared bright 
nebula, or unresolved infrared excess. The observed frequency of infrared cir- 
cumstellar emission (< 24%) and the mass loss rate we estimate for our sources 
shows that dusty mass loss can only account for part of the Cepheid mass loss 
discrepancy. Nevertheless, our direct evidence that mass loss is active during 
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the Cepheid phase is an important confirmation that these processes need to be 
included in evolutionary and pulsation models of these stars, and should be taken 
into account in the calibration of the Cepheid distance scale. 

Subject headings: Cepheids — infrared: stars — stars: mass loss 



Introduction 



Classical Cepheid variable stars are of key importance in the extragalactic distance 
scale and in studies of Galactic stellar and chemical evolution. But the very property that 
makes them useful as distance indicators — pulsation — also makes Cepheids difficult to model 
accurately. For several decades, a major problem in the understanding of Cepheids has been 
the discrepancy betwee n masses estimated from ste l lar evolution theor y and those estimated 
from pulsation theory (jChristyi Il968r IStobid Il969l : iFricke et al.l 119721 ) . While the problem 
was partly addressed by revision of radiative opacities and changes in evolutionary tracks 



(IMoskalik et al.l Il992l) . the discrepancy persists at the 10-20% level ( jCaputo et al.l 12005 



Keller fc Wood! l2006h . Measurement s of several Cephe i d mas ses are now available from 
binary systems ([Benedict et al.l 120071 : Evans et al.l 12008k 120091 ). but the comparison with 
evolutionary and pulsational masses has not yet produced a definitive result. 

The re are a number o f possible solutions to the 'Cepheid mass discrepancy' summa- 
rized by iBono et al.l (120061) . Recent studie s have come to differing conclusions on which 
might be responsi ble: ICaputo et al.l (120051 ) suggested that mass loss can account for the 
discrepancy, while iKellerl ( 120081 ) concluded that additional mixing in Cepheid main sequence 
progenitors was responsible. Mass loss remains one of the less-understood parameters in 
stellar evolutionary theory (jWillson et al.l 120081 : IVinkl 120081 ) . and its relevance to Cepheid 
masses is still not completely clear. In the last several years, several lines of investigation 
have suggested that mass loss may be important in Cepheids. These i nclude the discovery 



of ex t ended near-infrared emission from a number of Galactic Cepheids (IMerand et al.ll2006 



20071 : iKervella et al.ll2006l ). p ossibly from circumstell ar shells, as well as theoretical modeling 



Cepheids (INeilson et al.l 12009 



of pulsatio n-driven mass loss (Neilson fc Lesterll20081 ) and its application to Magellanic Cloud 



20101 ). Both the presence of shells around nearby Cepheids 



and the evidence for infrared excess in LMC Cepheids have implications for the stars' de- 
rived diameters and therefore for the distance scale. Systematic effects on distances are 
possible, if Cepheid mass-loss is substantial enough to affect spectral energy distributions, 
and if it work s diffe rently in nearby and distant galaxies (e.g., is a function of metallicity). 
Neilson et al.l ( 120101 ) concluded that mass loss in Cepheids "has an important effect on the 
structure of the IR P-L [infrared period-luminosity] relations" although the extent of the 
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effect is not well characterized. 



Attempt s to estimate Cephe i d mas s-loss rates through infrared observations began about 
25 years ago. iMcAlary &: Welch! (119861 ) used infrared excesses inferred from the IRA S Point 
Sourc e Catalog to estimate mass loss rates of M < 10~ 9 — 10~ 8 M yr -1 , while iDeasy 
(119881 ) combined ultra yiolet and infrared o bservations to estimate upper limits of M < 



1CT 10 — 10~ 7 M Q yr -1 . IWelch fc Durid (119881 ) used radio observations to place similar upper 
limits on the mass loss in ionized gas. In general these rates are too low to account for the 
mass discrepancy. Recently, the Spitzer Space Telescope has been used in a number of studies 
of Cepheids in the mid-infrared, including the per io d-luminosity relation (INgeow &: Kanbur 



ol Uepneids m tne mid-mlrared, including tne per io d-luminosity relation niNgeow &: Kanbur 
20081 : iFreedman et alJliooilMarengo et alJboiOah . iKervella et all (bopgh combined ground- 
based infrared imaging and interferometric measurements with Spitzer data (including data 
from the program described here) to detect compact circumstellar envelopes around both 
I Car and RS Pup; the latter was previously known to be associated with a reflection nebula. 

To investigate the poss ibility of infrared excesses due to mass loss, we have obtained 
Spitzer (jWerner et al.ll2004j ) observations of a sample of 29 Classical Cepheids. Our results 
on the mid-infrared perio d-luminosity and period-color relations for Cepheids are presented 



Marengo et al.l ( 



2010al. Paper I); a detailed analysis of extended emission near 5 Cep is 



in 

presented by iMarengo et al.l (j2010bl ). In Paper I we determined that there was no evidence 
from IRAC colors for warm circumstellar dust. The present paper investigates the pres- 
ence of cooler dust around this sample of Cepheids via a search for extended mid-infrared 
emission. An important consideration is that extended emission may be the result of recent 
star-forming activity near a Cepheid rather than mass loss from the Cepheid itself. However, 
distinguishing those two different causes is not straightforward, so a search for extended emis- 
sion places an upper limit on the amount of mass loss. Except whe re mentioned otherw ise, 
in all of the following analysis we use the stellar distances given by lFouque et al.l (120071 ). 



2. Observational data 



The sample selection and first-epoch observations are described in detail by lMarengo et al 



(j2010al ). Briefly, 29 Ga lactic Cepheid va riables and 3 non-variable su pergiants were ob- 
served with both IRAC jFazio et al.ll2004h and MIPS faieke et all I2004T ) as part of Spitzer 
Space Telescop e Gener al Observer, program ID 30666. The IRAC observations described by 
Marengo et al.l (j2010al ) were primarily made in subarray mode. The 6 brightest stars were 
observed in full-array mode so that fluxes could be measured by fitting the outer part of the 
point spread function (PSF); 4 stars were observed in both full-array and subarray modes. 
The MIPS observations used the Photometry Astronomical Observing Template with the 
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shortest possible observation time (48.2 s) at 24 fim and varying times (37.7-1300.2 s) at 
70 fim. Three stars (S Nor, V Cen, U Sgr) were observed in two adjacent fields of view at 
70 (Ltm to better cover the su rrounding clu sters. While our target SZ Tau (see §3.5p may also 
be in the cluster NGC 1647 (ITurnerlll992[ ). after inspecting its 24 /zm image, we changed the 
planned 70 fim observation to use only one field of view centered on the star. This allowed 
us to increase the observation time and, as a bonus, also provided a second epoch of 24 /im 
observations. We performed simple aperture photometry on post-BCD mosaics primarily 
from the pipeline version S16.1 output. For the three cluster stars, the 70 jum data were 
remosaiced using MOPEX to combine the data for the two fields of view. For two stars 
(RS Pup, GH Lup) the 70 jum mosaics were remade to remove negative sidelobes. 



2.1. New IRAC data and reduction 



To follow-up earlier results, we obtained some additional Spitzer data in Cycle 5, on 
the dates 2008 October 5-6. All stars were observed in full-array mode with IRAC as part 
of Guaranteed Time program 50350 (PI G. Fazio). This increased the field of view for 
stars which had previously been observed only in subarray mode, and also provided deeper 
observations (2-second frames at each of 36 dither positions) for all stars. Each target fainter 
than K = 2 was observed with the IRAC full frame 12 second High Dynamic Range (HDR) 
mode, using a 12 point Reuleaux small scale dither pattern for a total integration time of 
144 sec in each band. HDR 12 second frames were used to better reveal the area closer to the 
saturated cores of the targets, which also helped the PSF subtraction process. For the six 
brightest sources, in order to avoid excess saturation, targets were observed with 2 second 
frames using IRAC full frame non-HDR mode and the 36 point Reuleaux small scale dither 
pattern, for a total integration time of 72 s in each band. The Reuleaux dither pattern was 
used to achieve the required total exposure time while obtaining the best possible spatial 
sampling of the targets in anticipation of the PSF subtraction process. 

The basic calibrated dat a (BCD) for all obse rvations were reduced using the custom 
post-BCD software IRACproc (jSchuster et al.ll2006l ). which generated mosaics for each source 
with the final exposure depth and outliers removed. The final images had a pixel scale of 
0.8627 arcsec pix -1 (l/v2 the original IRAC pixel scale) which provides a better sampling 
of the sources' PSF for a more accurate PSF subtraction. The PSF subtraction was also per- 
formed usi ng a routine part of IRACproc, which has been ext ensively tested in other Spitzer 
programs (IMarengo et al.l 120061 . l2010al ; iLuhman et al.l 120071 ) to search for faint structures 
and point sources around bright stars. 
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2.2. MIPS photometry 



Although the MIPS observations used in this paper are the same as those in lMarengo et al 



(l2010af ). a major purpose of the present work is to examine extended emission around a sub- 
set of the stars in our sample. Accordingly, we have performed point source fitting and 
subtraction on the 24 /im images of some targets to separate the contributions of the star 
and extended emission to the total flux density. The data were reprocessed with the MIPS 



instrument team Data Analysis Tool using standard procedures. The PSF used 



or th e 24 /im 



subtraction was built from the observations of the A4 V star, T3 Eri dSu et al.l 120081). which 
matches well with the theoretical PSF computed from STinyTim Program ( )Kristll2006l ) after 
proper smoothing. The subtraction was done by iteratively shifting the PSF to match the 
centroid of the target in sub-pixel accuracy, adjusting its scaling, and subtracting it from 
the target image until the residual inside r < 12" (~ to the first bright Airy ring of the PSF) 
was minimized. The residual flux is generally less than 1% of the total flux of a typical point 
source on a clean background. 

To estimate the non-point source emission near stars showing extended emission in the 
MIPS bands, we used the results of point source fitting and subtraction. We performed aper- 
ture photometry of the point and extended emission on images which had the field (but not 
the target) stars subtracted, using the largest apertures which could be contained in the im- 
age. From this we subtracted the flux measured by fitting and subtracting the target stars to 
derive the total non-point source emission. This procedure depends on accurate background 
subtraction for the large-aperture photometry, particularly difficult in the 70 /im images, 
and the resulting uncertainties are large. An alternative approach which can cons t rain th e 



physical conditions within regions of extended emission is used by iMarengo et al.l ( I2010bl ): 
measuring the surface brightness of the extended emission in small areas. Because the stars 
are heavily saturated in the IRAC bands, this is the only method available to investigate 
the extended emission visible in the IRAC images. In both IRAC and MIPS images we 
measured the surface brightness in several small regions near stars with extended emission, 
and subtracted backgrounds measured in regions located off the extended emission. The 
regions used are shown in the images of the relevant stars; the measured surface brightnesses 
are reported in Table [TJ 



3. Analysis 

Infrared colors were used in Paper I to demonstrate a lack of warm circumstellar dust 
in our sample of Cepheids. However, dust mass lost from Cepheids and then transported to 
larger distances would be expected to be cooler, and possibly visible as extended infrared 
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emission. One examp le of this phenomenon is the infrared nebula visible around 5 Cep 



(jMarengo et al.ll2010bf ): here we complete the survey of our Cepheid and comparison star 
sample. None of the stars showed any extended emission in the 3.6 fim or 4.5 fim bands, 
so those bands are not discussed further here. Distinguishing between extended emission 
associated with a specific star and that due to nearby Galactic cirrus and/or star formation, 
especially at low Galactic latitudes, is not always straightforward. We expect that emission 
associated with an individual star should be centered on or near the star, with its surface 
brightness decreasing with increasing distance from the star. Unrelated infrared emission 
will not in general have these characteristics, although of course it could hide or confuse 
real low level emission. For brevity, we refer to unrelated infrared emission as 'cirrus' below, 
although we recognize that it may not be true Galactic cirrus emission. 

Of our sample of 29 Cepheids, 16 showed clear signs of cirrus in the IRAC 8 /im images: 
AQ Pup, BF Oph, DT Cyg, FF Aql, GH Lup, RS Pup, S Mus, S Nor, SW Vel, T Mon, 
U Car, U Sgr, V636 Sco, V Cen, VY Car, and X Cyg. Cirrus is also visible in the MIPS 
images of GH Lup, V Cen, VY Car and U SgrQ As expected, the presence of cirrus is 
related to Galactic latitude: of the Cepheids in our sample with |fe| < 10°, 15 of 23 appear 
to be surrounded by cirrus emission. The remaining star with apparent unrelated emission 
is DT Cyg, at b = —10.8°. The list of cirrus sources include three of the four stars we 
believed as cl uster members (S Nor, V Cen, U Sgr); this classification is confirmed by 



Gieren et al.l (119971 ). who also found that VY Car is a member of the Car OB 2 association. 
The extended emission near these stars could be associated with either their clusters or the 
general interstellar medium. Figure [T] shows some examples of stars surrounded by cirrus. 

As a quantitative test for excess infrared emission near the stars, we computed the aver- 
age surface brightness around each star and compared it to a background measurement. The 
IRAC 8 fim images were used for this measurement as they have a larger uniform-exposure 
field of view than the MIPS images, and much better sensitivity to extended emission than 
the 5.8 /im images. The average surface brightnesses were measured in 1 arcminute radius 
apertures surrounding each star, with the PSF subtraction residuals and bright star artifacts 
masked. The background measurement was made in an annulus with inner and outer radii 
of 74 and 98 arcseconds. Point sources near the target stars, imperfect PSF subtraction, 
and unrelated cirrus emission all complicate the comparison between low-surface-brightness 
emission near the stars and in their vicinities. To be conservative, we consider as having 
possible extended infrared emission only those stars for which the average surface bright- 
ness in the measurement aperture exceeds that in the background aperture by 20 times the 



1 Of our three non- variable comparison stars, HD 1822296 and HD 183864 were observed in subarray 
mode only; ip And was observed in full-array mode and shows no evidence for cirrus. 
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quadrature sum of the surface brightness standard deviations. Since our measurement errors 
are are influenced by the factors listed above and are unlikely to be Gaussian, we do not 
consider these detections to have true 20a significance — the value chosen is intended as a 
starting point to guide the analysis. 

We found no evidence for extended infrared emission around two stars which have been 



previously found to have circumstellar envelopes: Polaris (IMerand et al.ll2006l ) and Y Oph 



jMerand et alil2007h fl Other stars with neither cirrus or extended infrared emission are BB 



Sgr, U Aql, V350 Sgr, W Sgr, (3 Dor, r\ Aql and ( Gem. Five stars showed both cirrus and 
extended emission centered on or near the star: GH Lup, RS Pup, S Mus, T Mon, and X 
Cyg. Three further stars had little detectable cirrus but did have extended infrared emission 
centered on the star: 5 Cep, £ Car, and SZ Tau. 



The extended emission around 5 Cep is discussed in detail by iMarengo et al.l (l2010bl ) 
and briefly summarized here. The emission consists of a large-scale (0.1 pc) arched structure 
visible at 70 jttm, with diffuse, filamentary emission at 5.8, 8.0 and 24 /im contained within 
the arch. The brightest region of this diffuse emission is between S Cep and its hot companion 
HD 213307. The symmetry axis of the 70 fim emission near S Cep is aligned with the star's 
velocity vector relative to the local ISM, suggesting a bow-shock at least partly swept-up 
from the ISM. The emission near the stars is hypothesized to be due to the interaction of 
their stellar winds. The mass loss rates required for the creation of these structures are 
estimated to be in the range 5 x 10~ 9 — 6 x 10~ 8 M yr _1 . 

The remaining stars with extended infrared emission are discussed individually below. 



3.1. GH Lup 

Images of GH Lup are shown in Figure [21 The most obvious component of extended 
emission is an arc-like structure to the east of the star, visible at 8, 24, and 70 /im. The 
structure is not centered on the star itself but roughly 20" to the south (2.2 x 10 4 AU at 
the star's distance) with a radius of approximately 32" (3.6 x 10 4 AU). Spatially coincident 
with the arc is a point source visible at 5.8 and 8.0 /^m only, which we can identify with the 
K s = 9.5 object 2MASS J15244147-5251282. Because the region is crowded, with copious 
diffuse emission, it is unclear whether the arc near GH Lup is related to the Cepheid, 
the nearby star, or neither. The stellar crowding also makes defining regions for surface 



2 The extended emission found by those authors was at very small spatial scales (a few stellar radii), well 
below our angular resolution. 
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photometry of the diffuse emission difficult; we chose the regions marked in Figure [2] to 
avoid stars in the 8 /im image and capture at least some of the extended emission. The 
resulting measurements are given in Table [U with the colors of the extended emission plotted 
in Figure [3J Using a modified blackbody spectrum to assign color temperatures from the 
S24/ S70 ratios yields 55 < T < 81 K, depending on the emissivity index j3 for the modified 
blackbody The Ss/SW ratios for GH Lup tend to be higher than those found for the other 
stars with extended emission. Together with the structure of the emission and the extensive 
nearby cirrus, this leads us to speculate that the emission near GH Lup is only moderately 
likely to be associated with the star itself (see §@|. 



3.2. 



Car 



Images of I Car are shown in Figure HI These data were also analyzed by iKervella et al. 



( 120091 ) . who failed to detect any nebulosity in the Spitzer images, but found evidence of 
small scale (~ 1") extended emission in near- and mid-IR interferometric data, attributed 
to hot circumstellar dust. Our PSF-subtracted IRAC images show low surface brightness 
emission, at 5.8 an d 8.0 nm, surroundin g the star with a radius of ~ 35", or 1.7 x 10 4 AU at a 
distance of 498 pc ( IKervella et al.ll2009l ). We do not detect any significant extended emission 
at 70 /im, and only a tentative excess at 24 /im, northeast of the nearly saturated core of the 
PSF-subtracted star. The 8.0 /im emission appears brighter northwest of the star, where an 
arc-like structure is detected at a distance of ~ 30" (1.5 x 10 4 AU). Following a suggestion by 
the referee, we computed the components of the star's Galact ic peculiar space motio n and 
their projection onto the sky, using the procedure described bv lMarengo et al.l (I2010br). The 



star's position and prop er motion were taken from the Hipparcos catalog ( iPerryman et al. 



19971) , the distance from IKervella et al.l ( 120091 ). and heliocentric radial velocity of 3.6 km s 



from Nardetto et al. rl2009h . We find a space velocity of 24.9 km s 1 along a position angle 
of 315° (east of north), with the direction shown as an arrow in the second panel of Figure 
HJ The vector is roughly aligned along the direction from the star to the arc-like structure, 
suggesting a possible connection between the two (e.g., a bow shock, reflection nebula, or 
light echo). Image artifacts from the PSF subtraction and the star make the connection 
difficult to secure, however. Artifacts also make it impossible to accurately measure the 
surface brightness of the extended emission, and thus estimate the mass loss rates that could 
be responsible for the observed nebulosity. Higher spatial-resolution imaging covering a 
wider area would clearly be useful in further study of £ Car. 
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3.3. RS Pup 



RS Pup is embedded in a well-known reflection nebula dWesterlund Il961t lHavlenlll972 ) 

whos e infrared emission has been noted by numerous authors (e.g., Gehrzl 19721 ; McAlary fc Welch 
19861 ). PSF-subtracted images of RS Pup are shown in Figure [5j The morphology of the 
extended infrared emission varies with wavelength but the brightest region is to the north- 
east of the star. Most of the emission is within about 80" of the star, a linear distance of 



0.77 pc at the 1.992 kpc distance of RS Pup (IKervella et al.ll2008l ). although the small field 
of view of the 70 /xm image makes this diffic ult to pin down . The MIPS data, but not the 



full-frame IRAC data, were also analyzed by iKervella et a 



(20 091), so he r e we present the 



first image of the resolved 8.0 /im emission around RS Pup. IKervella et al.l ( 120091 ) attributed 
the > 20 /im resolved emission, as well as the infrared excess in the IRAS data, to an ex- 
tended envelope with a radius of roughly 1 pc and a temperature of ~ 40 K. IKervella et al. 



( 120091 ) also detected resolved emission at much smaller angular distances to the star, and 
attributed this to a "warm and compact component" at a distance of 300-3000 AU and a 
temperature of a few hundred K. They proposed that the warm component comes from RS 
Pup's mass loss, and the cold componen t is interstellar medium gas compressed by the star's 
stellar wind. In iMarengo et al.l ( 12 01 Obi ) a similar origin for the extended infrared emission 
around 5 Cep was proposed. 

We add to the analysis of RS Pup's extended emission by examining the colors of the 
extended emission. Using the surface brightness measurements in Table (H Figure [3] plots 
the ratios of the extende d emissi on around RS Pup and S Mus, with the values derived for 
5 Cep by IMarengo et al.l (l2010bf ) for comparison. Using a modified blackbody spectrum to 
assign color temperatures from the S24/S70 ratios yields 50 < T < 80 K, depending on the 
emissivity index for the modified blackbody. These temperatures are a little lower than 
the values derived for 5 Cep; however, compared to 5 Cep, the RS Pup nebulosity is more 
distant from the central star. Like the 5 Cep nebulosity, the RS Pup extended emission 
has S&o/Su ~ 0.5, suggesting that some polycyclic aromatic hydrocarbon s (PAHs) must 



be pre sent (pure dust at 50 < T < 80 K would imply a much lower ratio; IMarengo et al. 



2010bl ). We can also compute an average dust temperature for the RS Pup nebulosity using 
the total flux densities of the extended emission. At 70 /im we used the largest aperture 
that will fit on the image (85" radius) to derive a total flux density of 9.7 Jy. At 24 /im we 
measured the flux density in a large aperture (72" radius) and subtracted the PSF-fitting flux 
of the central source to derive an extended emission flux density of (620 — 330 =) 290 mJy 
Using a modified blackbody spectrum to assign color temperatures from the ratio of the two 
bands yields 44 < T < 57 K, depending on . Thes e temperatures are a little higher than 
the estimate of 40 K given by IKervella et al.l ( 120091 ) , who extrapolated the 70 /im surface 
brightness profile to derive a larger total flux density. 
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The dust mass in the RS Pup nebula is best estimated with far-infrared photometry, 
since as Figure [5] shows, the stellar emission is insignificant at wavelengths of 70 /im and 
beyond. Because our 70 /im MIPS images have a rather small field of view, the difficulty 
of accurate back ground subtraction m akes any attempt to estimate the total flux density 
(including that of iKervella et al .1120091 ) problematic. However, a rough estimate of total dust 
mass is possible and provide s a co mparison point for previous results. We use the formula 
given in Eq. 2 of Evans et al.l (120031 ). assuming that all of the 70 /im emission is from optically 
thin dust with T = 40 — 80 K and emissivity n v = 3Q„/pa in cm 2 g~ x . For comparison with 
McAlary fc Welch! (119861 ). we use k = 19.0, scaled from the 100 /im value the a = 0.1 /tm 



'dirty silicate' dust model in I Jones fc Merrill! (119761 ) via a Q u oc A~ 2 wavelength dependence. 
The derived total dust masses for RS Pup are in the range 1 x 10~ 3 — 1.4 x 10 -2 M Q . (The 
range of values here accounts for the change in the Planck function, but not any temperature 
dependence of the dust emissivity, over the given temperature range.) The uncertainty 
due to temperature is compounded by the uncertainty in dust emissi vity: for example , 
using the wavelength-dependent emissivity of amorphous fayalite given by lEvans et al.l (120031 ) 
would result in a value of k nearly 15 times larger (and thus masses 15 times smaller) 
than given above . The 100 /im flux density from the IRAS Small Scale Structure Catalog 
(jHelou &: Walker! Il988l ) should provide a better estimate of the t otal flux density of the 
RS Pup nebula. Using this measurement, iMcAlary &: Welch! (119861 ) estimated a total dust 
mass of 0.028 M (for their adopted distance of 1.778 kpc), or 0.035 M for the 1.992 
kpc distance we use. The mass derived from the 70 /im flux density is the same order of 
magnitude, likely smaller because our observations do not capture the full extent of the RS 
Pup nebula. 



As IMcAlary fc Welch! (119861 ) point out, the large size of the RS Pup nebula makes 
estimating a mass-loss rate from the nebular m ass "extremely difficult". With a gas-to- 
dust ratio of 100, the IMcAlary fc Welch! (119861 ) dust mass converts to a total mass of 
3 M . An argument in favor of the idea that some of the nebula is swept-up interstellar 
medium (ISM) results from a simple comp arison of the nebula mass with the helium-burning 
ti mescale of RS Pup ( Kervella et al.l 120091 ) . Taking the mass estimate for RS Pup provided 
bv lCaputo et al.l J20051 L the helium burning time inside the instability strip is of the order of 
1.2 Myr. If all of the nebula consists of mass lost by RS Pup during its helium-burning phase, 
the implied mass- loss rate w ould be > 10~ 6 M^yr -1 . This would imply a total mass loss at 



least of the order of 10%. As IKervella et al.l (120091 ) note, this would be suf ficient to solve the 
pulsat ional/evolutionary mass difference for RS Pup. However, as noted by lNeilson &: Lester 
(120081 ). the observational mass loss rate for RS Pup disagrees with their theoretical predic- 
tions (1.6 — 6.5 x 10 -9 M yr _1 ) by several orders of magnit ude. These predictions do appear 
to fit the observational data for Magellanic Cloud Cepheids (INeilson et al.ll2009l . l2010l ). which 
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makes it more believable that at least some of the mass in the RS Pup nebula did not orig- 
inally belong to the star. The observed S 8 , /S 2 4, ratio is too high for pure dust emission, 
but lower than typical for the typical PAH content of the ISM. Since Cepheid atmospheres 
are oxygen-rich and not expected to form PAHs, this suggests that both swept-up ISM (or 
pre-existing material containing PAHs) and stellar wind must comprise the nebula. 



3.4. S Mus 



Images of S Mus are shown in Figure There are slightly extended emission features 
to both the north and south of the star visible at both 8 and 24 /zm, with radial extents 
35-45" (2.9 - 3.7 x 10 4 AU at a distance of 820 pc). There is also extended emission visible 
at 70 /im, to a radius of about 55" (4.5 x 10 4 AU), although the morphology is different at 
that wavelength. The surface brightness of this emission is faint compared to that around 
RS Pup, and the orientation of the PSF subtraction residuals makes computing Sg.o/SW and 
S24/ S 70 at the same location difficult. The S Mus point shown in Figure [3] is for the box 
marked '3' in the IRAC images and '2' in the MIPS images. Although uncertain, the colors 
imply an even lower PAH content that in RS Pup and 5 Cep, and a som ewhat higher color 
temperature (70-120 K). This star has a high-temperature companion (jEvans et al.l 120061 ) 
and the extended emission could be associated with that star rather than S Mus itself, or 
the compa nion could be respons ible for heating the circumstellar material (as hypothesized 
for 5 Cep; IMarengo etaPhoiObh . 



The photometry of the extended emission around S Mus can be used to estimate the 
mass lost by the star. Large-aperture (35" radius) photometry of the 70 jum emission yields 
a total flux density of 96 ± 12 mJy; subtracting the point-source flux density (measured in a 
16" radius aperture as for Paper I) of 39 ±11 mJy yields a nebular flux density of 57±16 mJy. 
We use the same dust mass estimation procedure as for RS Pup, a temperatures of 70-120 
K, and a stellar distance of 0.88 kpc to derive an implied dust mass 4 — 17 x 10~ 7 M . With 
a gas-to-dust ratio of 100, the tot al implied mass would b e approximately 4 — 17 x 10~ 5 M . 
If we instead use k = 56 as in IMarengo et al.l ( )2010bl ). the total mass agrees quite well 
with that derived 8 Cep. Given the photometric uncertainties, the agreement may well be 
coincidental; however the order of magnitude comparison is intriguing. A mass loss timescale 
can be estimated by dividing the distance to the emission (35" is 0.15 pc at the distance 
of S Mus) by the escape velocity (100 km s _1 ) to yield t 1500 yr, so the resulting mass 
loss rate would be 3 — 10 x 10~ 8 M Q yr^ 1 . As for S C ep, this estimate i s con sistent with 



the predictions of iDeasyl ( 119881 ) but well above those of iNeilson &: Lesterl ( 120081 ). There are 
two reasons why our value for the mass loss rate is likely to be an upper limit, and therefore 
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larger than the iNeilson fc Lesterl (120081 ) predictions. (1) We assumed that all of the dust 
responsible for the extended emission comes from wind material as opposed to having been 
swept-up from the ISM; as for RS Pup, the PAH content implies that ISM sweeping is likely 
relevant. (2) Using the escape velocity to estimate the nebula's age results in the minimum 
possible age, and hence a maximum for the mass loss rate. 

We can estimate the total mass lost for S Mus through its mass and crossing time for 
the instability strip, as above for RS Pup. We es t imate d pulsational and evolutionary masses 
through the rela tion provided by ICaputo et al.l ( 120051 . Tables 2-5); using the n ear-infrared 



photo metry from lLaney Sz Stobie (|l994[ transformed to the 2MASS syst 



20071 ) and optical phot ometry from 



em using 



Koen et al. 



Pedicelli et al.l (120101 ) . correcte d for extinctio n usin g 



the reddening given by iFernie et al.l ( 119951 ) and the reddening law of ICardelli et al.l ( 119891 ). 
We used an intrinsic distance based on the near-infrared photometry and period- Wesenheit 
relations provided by Persson et al. ( 2004j) . This procedure results in a mean mass 7.0 ± 



1.5 M Q , and according to lBono et al.l (120001 ) the corresponding crossing time of the instability 



strip is of the order of 0.14 Myr. Given the mass loss rate above, this results i n a small total 
mass loss of 0.006 M . Using the recent estimate of the radius of S Mus by iGroenewegen 
(120081 ) and the period- mass-radius relation given by lBono et al.l (120011 ). the pulsation mass is 
lower, ~ 5 M , and the instability strip crossing time much longer, of order 1 Myr; however 
the total mass lost would still be only 0.04 M , too small to significantly affect the mass 
discrepancy. 



3.5. SZ Tau 



SZ Tau was the only Cepheid found in Paper I to have a statistically significant [24] — [70] 
color excess. Images of this star are shown in Figure [71 There is an extended object about 
25" x 8" clearly visible in the PSF-subtracted 8.0 /zm image, about 7" to the north-east 
of the star. At the 512 pc distance to the star, 25" x 8" corresponds to linear scales of 
(1.3 x 0.4) x 10 4 AU. The object is not clearly visible in the 5.8 /im band, due to PSF 
subtraction residuals, however it is visible at 24 /im and appears to be the source of most 
of the 70 fim. luminosity. The object is resolved at 24 /xm and has a similar appearance in 
two images in that band taken a bout a year apart, suggesting that it is not an instrumental 
artifact or a solar system object. iNoriega-Crespo et al.l (119971 ) detected a "ridge or filament" 
of similar shape but much larger angular size in IRAS imaging of Betelgeuse; they su ggested 
that this structure is emission from the interstellar medium, unrelated to the star. Turner 



(119921 ) concluded that SZ Tau is a first-overtone pulsator and noted several other unusual 
characteristics of this star (e.g., period changes), but it is unclear how these features could 
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be related to the presence of extended emission. 



The elongated shape of the emission near SZ Tau suggests a background galaxy. The 
infrared colors should yield further information on the object's nature, although photometry 
of the SZ Tau extended emission is complicated by the star itself. At 70 /im, the extended 
emission appears to be the dominant source of emission. When we measure flux density in 
a 35" radius aperture centered on the 70 /mi peak and a 39-65" background annulus, the 
resulting flux density (which may include some contribution from the star) is 74 mjyjj We 
measured the flux density of the SZ Tau extended emission at 24 /im by doing a large-aperture 
(35" radius) measurement on the extended emission+star and subtracting the PSF-fitting 
star flux, to give an extended emission flux density of (175 — 144 =) 31 mJy. Measuring 
the 8 /im flux density is the most difficult since there is significant contamination from the 
PSF subtraction residuals. Measuring the emission on an image with the residuals masked 
in 35" radius aperture with 35-43" background annulus gives a flux density of 10 mJy; 
applying the IRAC extended source photometric correction changes this to 7.4 mJy. Different 
approaches to residual masking and photometric aperture change the resulting flux density 
by ~ 20 — 30%, so the value is highly uncertain. This is hardly surprising since the extended 
emission is more than 100 times fainter than the star at 8 /im. Consistent with this result is 
the classification of the star's mid- infrared spectrum as normal, with no dust, based on data 
from the PHT-S instrument on the Infrared Space Observatory (jfiodge et al.l I2004T ) . 



We can ask whether the the galaxy-like appearance of the extended emission near SZ 
Tau is consistent with its mid-infrared fl ux densities and co lors. The object is rather bright 
for a galaxy: number counts at 24 /im Jshupe et al.l l2008h show that there are only 3.5-4 
galaxies per square degree brighter than about 15 mJy. Taken together, all of our 24 /im 
Cepheid images cover about 0.5 square degree, so we might expect 1 or 2 such galaxies 
in our survey. Figure [8] compares the colors of the SZ Tau extended emission with some 
comparison samples. The color /s.n//24 = 0-24 is consistent with the c olors of both nearby 
and distant galaxies (IMunoz-Mateos et al.l l2009i ; lHainline et al.l 120091 ) . However, the flux 
ratio /24//70 = 0.42, while inconsistent with an unobscured stellar source, is ra ther high for 



a gal axy detected at 70 /im; such galaxies typically have 0.05 < /24//70 ^ 0.11 (IFrayer et al. 



20061 ). The flux ratios and colors are m ore consistent wit h those of candidate young stellar 
objects detected in the Cepheus Flare (IKirk et al.l l2009l ) . although YSOs do not typically 
have lenticular morphologies. Planetary nebulae can have asymmetric shapes, but the SZ 
Tau emission has quite different colors (see Figure [8]). The best candidates to explain the 



3 The flux density for this star reported in lMarengo et al' ( 2010a ) is slightly different because that value 
was measured centered on the 24 ^tm star position. Also, there is a bright far-infrared source about 80" to 
the northwest of the star's position which may also contaminate the photometry. 
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extended emission near SZ Tau are either an unusually-shaped Galactic object (possibly a 
YSO) or a background galaxy with very weak cold dust emission. This mysterious object 
warrants further investigation. A more definitive conclusion on its nature could be reached 
with higher spatial resolution infrared observations, in which the flux from the star and 
extended emission were more clearly separated. 



3.6. T Mon and X Cyg 



These stars show extended emission in the IRAC 8.0 /im band with only tentative 
detections in the MIPS bands; images can be seen in Figures M and [TDJ In the image of T 
Mon, there appears to be cirrus across the entire field of view. However, there is brighter 
nebulosity closer to the star, specifically to the north and south in an hourglass-like shape 
extending to about 45" radius (6.2 x 10 4 AU, or 0.3 pc, at a distance of 1389 pc). The 
emission is centered around the star, suggesting that it may be locally illuminated. The 
extended emission near X Cyg has the second-lowest significance (after GH Lup) of the 
stars with emission above background levels. It is brightest to the northwest of the star and 
concentrated within a radius of 30" (3.5 x 10 4 AU, or 0.16 pc, at a distance of 1163 kpc). 
For both stars, the 24 /im images show hints of enhanced brightness at similar positions to 
the 8 /zm emission, but these positions are too close to the PSF subtraction residuals to 
make definitive conclusions, or to measure a surface brightness. There is no corresponding 
detectable emission at 70 /im. With both stars having extended emission observed at only 
one wavelength, the detections must be regarded as tentative. Observations with the James 
Webb Space Telescope, particularly its coronographs, would provide an avenue to confirm 
or refute them. If the detections are confirm ed, the properties of the se stars may provide 
a challenge in relating them to any mass loss. iNeilson fc Lester! (120081 ) predicted that mass 
loss rates, while showing large scatter with per iod, are generally l arger for long-period stars. 
T Mon is a long period star with a companion (jEvans et al.lll9991). wh ile X Cyg has a period 
of 16.4 d and no known companion hotter than type A5 (jEvanslll992l ). However, a low-mass 
companion with an orbit in the plane of the sky is not ruled out for X Cyg and might be 
important in forming and/or shaping an extended nebula. 



4. Discussion and Conclusions 



Spitzer mid-infrared observations of C epheids have produce d no strong evidence for 
large amounts of warm circumstellar dust (jMarengo et all l2010al ). Here we use the same 
observations to place limits on cooler dust. Of our sample of 29 Cepheids, two stars (RS Pup 
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and 5 Cep) have well-detected extended infrared emission, two (GH Lup and S Mus) have 
possible extended emission at 24 and 70 /zm, and a further three (£ Car, T Mon and X Cyg) 
show tentative evidence for extended infrared emission at shorter infrared wavelengths only. 

Is there anything particularly special about the stars in which we do have (tentative) 
evidence for mass loss? Our observations of the nebulosity around S Cep and S Mus sug- 
gest that the detectability of a dustless wind might be enhanced when interstellar matter 
around the Cepheid is swept up into a bow shock, and a hot companion heats the material. 
Compared to the emission around RS Pup and GH Lup, the S24/SV0 ratio (a proxy for dust 
temperature) is higher near S Mus and 5 Cep. Both of these stars have known hot compan- 
ions, while RS Pup and GH Lup do not (although low-temperature companions cannot be 
excluded using current data). A counter example might be provided by T Mon, which also 
has a hot companion, but does not have detectable extended emission at 24 or 70 /im or an 
infrared excess; however, T Mon's companion is cooler (type AO) than that of the other two 
stars (B3 and B7-8). 

The surface brightness ratios Ss. 0/^24 are a rather crude tool for deriving detailed phys- 
ical conditions, but can give some insight into the PAH content of the ISM near these stars. 
The low PAH content indicated by the Ss.o/S^ ra ti° s f° r RS Pup, S Mus and 5 Cep — much 
lower than in the NGC 7023 reflection nebula, or the ISM in general — suggests that the ISM 
around these stars has been diluted to some extent by the stellar winds. The higher PAH 
content near GH Lup would imply less dilution, suggesting that the emission there may be 
unrelated to the star. We caution that, while the amount of dilution might be indicated 
by the Sg. o/SW ratio, high spatial-resolution spectroscopy in which the stellar and extended 
emission are separated would provide much greater physical insight into the conditions near 
these stars. 

The IRAC-only extended emission near £ Car, T Mon, and X Cyg also suggests that 
these stars are somehow unusual. However, we have been unable to find any clear distinction 
between these stars and the three which also show extended MIPS emission, or between all 
seven stars and our full Cepheid sample. Like RS Pup, I Car and T Mon have long periods 
(> 27 d) and low temperatures; GH Lup, S Mus and S Cep have periods < 10 d, and X Cyg 
i s intermediate with a period of 16.4 d. Six of the seven stars have s uper-solar metallicity 



(IRomaniello et al.ll2008 



Andrievsky et al.ll2005uKovtyukh et al.ll2005l ) but the difference in 



mean metallicity betwe en these stars and the other stars in the sample is barely statistically 



significant (p = 0.025). lEvand ( 119921 ) placed strong limits on the presence of hot companions 
in Cepheids: of the stars in our sample, 3 of the 7 with possible extended emission have such 
companions (S Mus, 5 Cep, T Mon), as do 5 of the 22 stars without extended emission (an 
additional 3/23 are radial- velocity binaries). With such small numbers it is not yet possible 



- 16 - 



to say whether the presence of hot companions enhances the detection of extended infrared 
emission. 

RS Pup is clearly different from all of the other Cepheids in our sample in terms of the 
extent and brightness of its infrared emission. One possible explanation is that this star has 
swept up more material, either by having been on the instability strip for longer (unlikely, 
since long-period stars evolve faster) or by being located in a region of higher density in 
the interstellar medium. An alternative explanation is suggested by visual examination of 
a Schmidt plate taken by one of us (DW) in the region of RS Pup. The image reveals a 
number of near-identical dust clumps, with similar angular sizes, in the vicinity of the RS 
Pup nebula. If these clumps are at the same physical distance and are physically related, but 
only one of them has a Cepheid in it, it could be argued that they are all likely pre-existing 
nebulae — that is, not the result of Cepheid mass loss. The counter-argument would be that 
the ring-like s tructures observed in the surrounding nebula appear to be centered on RS Pup 



(jHavlenl 119721 ): this would be an extreme coincidence if the nebula were unrelated to the 



star. 

Before our Spitzer observations, only two Cepheids had unambiguous evidence for cir- 
cumstellar nebulae: i Car and RS Pup. We have now doubled that number with the addition 
of S Mus and 5 Cep; if extended emission around GH Lup, T Mon and X Cyg is confirmed, 
this would increase the sample size even more. A rough estimate of the fraction of Cepheids 
with substantial dusty mass loss can be made as 4 to 7 of 29, or ~ 14 - 24% The true 
fraction of Cepheids with mass loss could be below this limit if the observed nebulosity is 
related to the gas and dust from which the star formed. This is however unlikely, because 
even for small space velocities (~ 1 km s" 1 , an order of magnitude less than the observed 
velocity of 5 Cep and I Car), in 50 Myr (the time it takes for an intermediate mass star to 
become a Cepheid) these stars will have moved away from their birthplace by several tens 
of parsecs. This is much larger than the size of the observed nebulosity around the stars in 
our sample. 

Addressing the initial motivation for this study, it appears that mass loss, at least its 
dusty component visible in the infrared, can only partially account for the Ce pheid mass loss 
discrep ancy (e.g. up to l/5th of a 10% discrepancy in the case of 5 Cep, see iMarengo et al. 
2010b|). Our estimates of Cepheid mass loss (in the range 10 -8 to 10 -9 M yr - 1 ) are compara 



ble to previous studies based on the detection of infrared and UV excess (e.g. IWelch fc Duric 



4 The shells detected in the near- infrared around some Cepheids ( Merand et al. 20061 2007 : Kervella et al 



2006) arc not consistent with dusty mass loss: t hey are very close to the stars and thus well above the dust 



sublimation temperature ([Marengo et al.ll2010al ). 
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19881 ; iDeasyi Il988h . but the higher spatial resolution of our images helps to separate local 
emission from background contamination. The detection of shocked structures, as in the 
case of S Cep and possibly i Car, is a compelling indication of current mass loss in those 
stars. While the total mass lost over the Cepheid phase is small, the confirmation that these 
processes are active needs to be taken into account for a proper modeling of Cepheid evolu- 
tion. The presence of significant Cepheid winds can affect hydrodynamic models of Cepheid 
pulsations, and the way shock waves propagate in the Cepheid atmospheres. Finally, the 
observed level of circumstellar dust needs to be f urther explored , as it may have a significant 



effect on inferred Cepheid distances as shown by iNeilson et al.l ( 120101 ) 
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Table 1. Surface photometry of extended infrared emission 



Star /region 


5.8 /im 


8 /im 


24 


70 jum 


GH Lup 1 


0.07 ±0.01 


0.27 ±0.01 


0.16 ±0.02 


1.1 ±0.1 


GH Lup 2 


0.02 ±0.01 


0.36 ±0.01 


0.29 ± 0.02 


1.8 ±0.1 


GH Lup 3 


0.13 ±0.01 


0.49 ±0.01 


0.51 ±0.02 


3.5 ±0.2 


GH Lup 4 


0.10 ±0.01 


0.41 ±0.01 


0.46 ± 0.02 


2.6 ±0.1 


RS Pup 1 


0.24 ±0.03 


0.63 ± 0.02 


0.89 ± 0.39 


9.4 ±0.8 


RS Pup 2 


0.18 ±0.01 


0.56 ±0.01 


0.82 ± 0.30 


24 ± 1 


RS Pup 3 


0.27 ±0.01 


0.83 ± 0.03 


1.99 ±0.86 


40 ±1 


RS Pup 4 


0.33 ±0.01 


0.91 ±0.01 


3.10 ±0.72 


47 ± 1 


RS Pup 5 


0.23 ±0.01 


0.84 ±0.03 


3.14 ±2.08 


35 ±2 


S Mus l a 


0.05 ±0.01 


0.10 ±0.01 


1.3 ± 1.1 


2.0 ±0.2 


S Mus 2 


0.08 ±0.01 


0.09 ±0.01 


0.8 ±0.6 


1.9 ±0.1 


S Mus 3 


0.15 ±0.01 


0.16 ±0.01 


0.7 ±0.3 


0.9 ±0.1 


S Mus 4 


0.07 ±0.01 


0.11 ±0.01 







For S Mus, regions are not the same in IRAC and MIPS. 



Note. — All surface brightness values in MJy sr 
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Fig. 1. — Images of some Cepheids with nearby extended infrared emission unrelated to the 
star. Left to right: DT Cyg, FF Aql, V Cen, and VY Car. Top row: PSF-subtracted IRAC 
8.0 jum images, Bottom row: non-subtracted MIPS 24 /im images. All fields of view are 4 
arcminute squares. In these and all following images, north is oriented up and east left. The 
black and white bands in the IRAC images are artifacts from both the bright source and the 
PSF-subtraction process. 




Fig. 2. — PSF-subtracted images of GH Lup in the (left to right) 5.8 /im and 8.0 /im IRAC 
bands, and the 24 /im and 70 /im (not PSF-subtracted; the white X marks the position 
of the star) MIPS bands. The field of view is a 4 arcminute square. Numbered boxes are 
regions used for surface photometry (see §2.21) . 
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Fig. 3. — Left: Surface brightness ratios of extended emission near GH Lup, RS Pup and 
S Mus, with emission near 5 Cep (circles) and NGC 7023 (star) shown for comparison. 
Increasing values on the horizontal axis correspond roughly to increasing dust temperature, 
while increasing values on the vertical axis correspond to increasing PAH content. 




Fig. 4. — PSF-subtracted images of I Car in the (left to right) 5.8 fim and 8.0 /im IRAC 
bands, and the 24 /im and 70 fim (not PSF-subtracted) MIPS bands. The field of view is a 
2 arcminute square. The white arrow in the second panel shows the direction of the star's 
proper motion relative to the local ISM. 



-25 - 




Fig. 5. — PSF-subtracted images of the reflection nebula around RS Pup in the (left to 
right) 5.8 /im and 8.0 /im IRAC bands, and the 24 /zm and 70 /im MIPS bands (not PSF- 
subtracted; the white X marks the position of the star). The field of view is a 4 arcminute 
square. Numbered boxes are regions used for surface photometry (see §2.2p . 




Fig. 6. — PSF-subtracted images of S Mus in the (left to right) 5.8 fim and 8.0 /im IRAC 
bands, and the 24 /im and 70 /tm (not PSF-subtracted) MIPS bands; the white X marks the 
position of the star. The field of view is a 3 arcminute square. Numbered boxes are regions 
used for surface photometry (see §2.2p . 
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Fig. 7. — PSF-subtracted images of SZ Tau in the (left to right) 5.8 jum and 8.0 /im IRAC 
bands, and the 24 /xm and 70 /im (not PSF-subtracted) MIPS bands; the white X marks the 
position of the star. The field of view is a 2 arcminute square. Of particular interest is the 
elongated object to the northeast of the star. 
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Fig. 8.— Right: flux de nsity ratios for subm i 



2009) , nearby g alaxies ( jMunoz-Mateos et al. 



limet er-detected galaxies (lHainline et al. 
2009h. Cep heus young stellar objects 



(IKirk et al.l 120091 ). LMC planetary nebulae ( jHora et al.l 120081 ). and the extended emission 
near SZ Tau. 



Fig. 9.— PSF-subtracted images of T Mon, in the (left to right) 5.8 /im and 8.0 fxm IRAC 
bands, and the 24 /im and 70 /im (not PSF-subtracted) MIPS bands, All images are 2.5 
arcminutes square. 




Fig. 10. — PSF-subtracted images of X Cyg in the (left to right) 5.8 /im and 8.0 /jm IRAC 
bands, and the 24 /im and 70 /im (not PSF-subtracted) MIPS bands. The field of view is a 
2.5 arcminute square. 



